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I. INTRODUCTION 
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In this paper, we study the degeneracies among several cosmological parameters in detail and 
discuss their impacts on the determinations of these parameters from the current and future obser- 
vations. By combining the latest data sets, including type-la supernovae "Union2.1" compilation, 
. . WMAP seven-year data and the baryon acoustic oscillations from the SDSS Data Release Seven, we 

04 ■ perform a global analysis to determine the cosmological parameters, such as the equation of state of 

dark energy w, the curvature of the universe ilk, the total neutrino mass ^ m^, and the parameters 

■ associated with the power spectrum of primordial fluctuations {us, Oa and r). We pay particular 
attention on the degeneracies among these parameters and the influences on their constraints, by 
with or without including these degeneracies, respectively. We find that w and ilk or ^ m,^ are 

■ strongly correlated. Including the degeneracies will signiflcantly weaken the constraints. Further- 
\^ \ more, we study the capabilities of future observations and flnd these degeneracies can be broken 

. very well. Consequently, the constraints of cosmological parameters can be improved dramatically. 

u: 
^' 

C^ Measuring cosmological parameters with observational data plays a crucial role in establishing the modern cos- 
Q . mology. In the past few years, with great progresses on astronomical observations, such as type-la supernovae (SN) 
?H ■ [H, cosmic microwave background radiation (CMB) [2| and large scale structure (LSS) (Sl, a standard cosmological 
c/3 . model is build up. It is commonly believed that our universe has undergone an inflationary period at the beginning 
and is accelerating expansion at present. Dark energy, the mysterious source jiriving the present acceleration of our 
] universe, has been studied widely in the literature since its discovery in 1998 [J]. However, the nature of dark energy, 
^-H ■ encoded in its equation of state (EoS) parameter w, remains controversial. For the other acceleration in the very 
^ ] early universe, the mechanics of inflation can naturally explain the flatness, homogeneity and isotropy of our universe, 
t . Inflation stretches the primordial density fluctuations and seeds the presently observed CMB and LSS. Understanding 

■ the dynamics of inflation and the origin of late-time acceleration are big challenges in physics and astronomy. 

, With the accumulation of observational data from CMB, LSS and SN observations and the improvements of the data 
quality, there are a lot of efforts in the literature on constraining the cosmological parameters from various observations 
! 0) SQ) such as the EoS of dark energy models w, the curvature of the universe flk, the total neutrino masses 
^ m^, and those associated with the running of the spectral index and gravitational waves. However, using different 
input theoretical model can affect constraints of cosmological parameters, since some of cosmological parameters are 
correlated with each other. Including or ignoring these degeneracies will give signiflcantly different constraints on 
cosmological parameters. For example, the EoS of dark energy w correlates with the matter energy density Vim and 
the curvature through the geometrical distance relations. Due to the strong correlation between w and fifc, we 
obtain different constraints on dark energy EoS w in the framework of flat and non-flat universe, respectively 0- The 
, EoS of dark energy also anti-correlates with the total neutrino mass, namely, the larger neutrino mass can be mimicked 
by the more negative w for the same distance relations (8l-ITo|. In the framework of dynamical dark energy model, 
the constraint on ^ is signiflcantly weakened, when comparing with that obtained from the standard ACDM 
model. Furthermore, the tensor perturbation mode and the dark energy component are correlated, since they mostly 
affect the large scale (low multipoles) temperature power spectrum of CMB. Thus, including the tensor fluctuations 
in the analysis will change the constraint of w. Moreover, some other cosmological parameters, such as, the initial 
conditions, inflationary parameters, can also influence the numerical results through the degeneracies. 

At present, we have a lot of cosmological observations, such as several hundred SN samples, CMB temperature and 
polarization fluctuation, LSS surveys and so forth. It is possible for us to study the determination of cosmological 
parameters precisely, especially the parameters correlating with w. Thus, in this paper, we study the uncertainties of 
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the constraining cosmological parameters in detail, due to the correlations among them. By employing the Markov 
Chain Monte Carlo method, we combine the latest observational data, such as seven-year WMAP data (WMAP7), 
baryon acoustic oscillation (BAO) and the "Union2.1" compilation SN sample to constrain cosmological parameters 
in different input cosmological models, such as the non-flat universe, the massive neutrino mass models, and so on. 
We pay particular attention on the degeneracies among these parameters. Furthermore, we also generate the mock 
data for the future observations, such as PLANCK [ll| and Euclid [l^.[T3l| projects, as well as future SN observations 
[3] to study their capabilities on breaking these degeneracies. 

The structure of the paper is as follows: in section|lT]we describe the method and the current and future observational 
data sets we use; Section [IIII contains our main numerical results on the cosmological parameters, and the last section 
is the summary. 



II. METHOD AND DATA 



A. Method 



In order to get constraints in different cosmological models, we employ the modified version of the MCMC package 
CosmoMC 15] and perform the global fitting analysis with the current and future observational data. We parameterize 
the cosmology with the following parameters: 



(ujb, LOc, ^k, 6s, r. Wo, Wa,U, Us, As, as, r) 



(1) 



where cjb = Q.hh? and uJc = i^ch^, in which Q,b and ftc are the physical baryon and cold dark matter densities relative 
to the critical density, flk is the spatial curvature and satisfies flk -\- i^b + -\- ^dc = 1, 6s is the ratio (multiplied by 
100) of the sound horizon to the angular diameter distance at decoupling, t is the optical depth to re-ionization, 
is the dark matter neutrino fraction at present (l6j . namely. 



~ Pm ^ 93.14 eV ^Imh? 
The primordial scalar power spectrum Vs{k) is parameterized as [Tt} : 



(2) 



ln7',(A:) = lnA,(fc,o) + K(fcso)-l)ln( — j +y 



ln( A 

fcsO 



(3) 



where Ag is defined as the amplitude of initial power spectrum, measures the spectral index, Ois is the running of 
the scalar spectral index and r is the tensor to scalar ratio of the primordial spectrum. For the pivot scale we set 
kso = 0.05 Mpc~^ which is the default value of CosmoMC package, wo and Wa are the parameters of dark energy, and 
the EoS is parameterized as [ll| : 



Wdc(a) = Wo + Wa(l - a) 



(4) 



(noted as "CPL"), where a = 1/(1 + z) is the scale factor and it covers ACDM model by taking wq = — 1 and Wq = 0. 
As we know, for time evolving dark energy model, it is crucial to include dark energy perturbations in the global 
fitting strategy to constrain cosmological parameters (19I - I21I ]. In this paper we use the method provided in refs.|2ll[22l| 
to treat the dark energy perturbations consistently in the whole parameter space in the numerical calculations. 



B. Observational Data 



1. Current Data 



In our analysis, we consider the following cosmological probes: i) power spectra of CMB temperature and polar- 
ization anisotropics; ii) the baryon acoustic oscillation in the galaxy power spectra; iii) measurement of the current 
Hubble constant; iv) luminosity distances of type la supernovae. 

To incorporate the WMAP7 CMB temperature and polarization power spectra, we use the routines for computing 
the likelihood supplied by the WMAP team [l^. The WMAP7 polarization data are composed of TE/EE/BB power 
spectra on large scales (2 < ^ < 23) and TE power spectra on small scales (24 < I < 800), while the WMAP7 
temperature data includes the CMB anisotropics on scales 2 < £ < 1200. 
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Baryon Acoustic Oscillations provides an efficient method for measuring the expansion history by using features in 
the clustering of galaxies within large scale surveys as a ruler with which to measure the distance-redshift relation. 
The BAO aries because that the coupling of baryons and photons by Thomson scattering in the early universe allows 
acoustic oscillations in the plasma at early time. After photons decoupling, the acoustic waves can propagate a certain 
distance, which becomes a characteristic comoving scale and remains in the distribution of galaxies. It provides a 
particularly roburst quantity to measure j23 - l28j . BAO has been detected in the current galaxy redshift survey data 
from the SDSS DR7 [3']. It measures not only the angular diameter distance, Dyi(z), but also the expansion rate of 
the universe, H{z), which is powerful for studying dark energy [29j . Since the current BAO data are not accurate 
enough for extracting the information of Da{z) and H{z) separately [s^l, one can only determine an effective distance 
3^ : Dy{z) = [(1 + z)^ D\{z)cz / H {z)]^^^ . In this paper we use the BAO measurement given by the SDSS-II survey 

In our analysis, we add a Gaussian prior on the current Hubble constant given by ref. (ssj : Hq — 74.2 ± 3.6 
kms~^Mpc~^ (68% C.L.). The quoted error includes both statistical and systematic errors. This measurement 
of Hq is obtained from the magnitude- redshift relation of 240 low-z Type la supernovae at z < 0.1 by the Near 
Infrared Camera and Multi-Object Spectrometer (NICMOS) Camera 2 of the Hubble Space Telescope (HST). This is 
a significant improvement over the previous prior, Hq — 72 ±8 kms""'^ Mpc""'^, which is from the Hubble Key project 
final result. In addition, we impose a weak top-hat prior on the Hubble parameter: Hq G [40, 100] kms^^ Mpc^^. 

Finally, we include data from Type la supernovae, which consists of luminosity distance measurements as a function 
of redshift, Dl{z). In this paper we use the latest SN data sets from the Supernova Cosmology Project, "Union 
Compilation 2.1", which consists of 580 samples and spans the redshift range < z < 1.55 [l|. This data set also 
provides the covariance matrix of data with and without systematic errors. In order to be conservative, we use the 
covariance matrix with systematic errors. When calculating the likelihood from SN, we marginalize over the absolute 
magnitude M, which is a nuisance parameter, as done in refs. [sj]. 



2. Future Measurements 



Since the current observations can not give the conclusive conclusion, we also use the simulated data from future 
observations of Planck and Euclid, as well as the SN sample. 

For the CMB simulation we consider a simple full-sky (/sky — 1) simulation at Planck-like sensitivity. We neglect 
foregrounds and assume the isotropic noise with variance Nj"^ = = Nf^ /2 = 3 x and a symmetric 

Gaussian beam of 7 arcminutes full-width half- maximum (FWHM). We use the simulated Ci up to £ = 2500 for 
temperature ^ and £ — 1500 for polarization. 

For the future LSS survey, we simply consider the Euclid-like survey which will measure ~ 10^ galaxies with the 
redshifts z < 2. In the measurements of large scale matter power spectrum of galaxies there are generally two 
statistical errors: sample variance and shot noise. The uncertainty due to statistical effects, averaged over a radial 
bin Afc in Fourier space, is (ssj : 

^ ^ 4^^ + ^) ■ 

The initial factor of 2 is due to the real property of the density field, Vsmvoy is the survey volume, and n is the mean 
galaxy density. In our simulations for simplicity and to be conservative, we use only the linear matter power spectrum 
up to k max — 0.1/iMpc-\ 

For the future supernovae observation, we simulate 4000 SN samples distributed in 17 bins from z — 0.1 to z = 1.7. 
We also include 300 low-z SN from the Nearby Supernova Factory [s^ to improve the constraints. For all the 
supernovae samples, we assume the variance u = 0.15 in each redshift bin. 



We also use a more conservative choice £max = 1500 to simulate the mock temperature power spectrum and find that the obtained 
constraints on cosmological parameters are almost unchanged. 
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TABLE 1. la constraints on some cosmological parameters from Union2.1+WMAP7+BAO+HST with difTerent input 

theoretical models. 



models 


Wo 
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-0.94 ±0.18 


-0.40 ± 0.78 
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-0.93 ±0.19 
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FIG. 1: 2-dimensional cross correlation between wq and Wa in different input cosmological models: the black solid lines are 
given by the standard CPL model, the red dashed, blue dash-dot and green dotted lines are given by CPL © flk , CPL © , 
and CPL © r models, respectively. 



III. NUMERICAL RESULTS 



A. Current Data 



In this section, we present the numerical results from current observational data sets with different input theoretical 
models. Here, we mainly focus on the EoS of dark energy, the primordial power spectrum index, and the total neutrino 
mass. 



1. Dark Energy 



To study the correlations between the EoS of dark energy and other cosmological parameters, such as the curvature 
of universe, the total neutrino mass, and tensor perturbation mode, we perform several global analyses with different 
input cosmological models, for example, the flat universe with massless neutrino (noted as "CPL"), the non-flat 
universe with massless neutrino (noted as "CPL © fife"), the flat universe with massive neutrino model (noted as 
"CPL © nil,") and the flat universe with tensor perturbation (noted as "CPL (Br"). We summarize the constraints 
on Wo and Wa in table I. 

In figure [1] we plot the 2-dimensional constraints of wq and Wa with different input cosmological models. The 
upper left plot is obtained from the standard "CPL" model. The 68% C.L. constraints are wa = —0.94 ± 0.18 and 
Wa = —0.40 ± 0.78, which is consistent with previous works [2]. However, if we choose a different input cosmological 
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model, the constraints will change obviously. In the upper right panel of figure [1] we compare the constraints on wq 
and Wa from the flat (black solid lines) and non-flat (red dashed lines) universe cases, respectively. Apparently, the 
model with non-zero ftk give weaker constraints on dark energy EoS parameters and the median values are shifted, 
namely, the constraints are wq = —0.93 ± 0.19 and Wa — — 0.47ig'g7 at the 68% confidence level. We know that the 
dark energy and the curvature are correlated via the distance relation. The effect on the distance of dark energy EoS 
can be mimicked by the non-zero flk- Therefore, the degeneracy between them should be very strong. In our analysis, 
we find that the correlation coefficients are —0.32 and 0.58 between flk and wq or Wa, respectively. 

Then, we consider the effect of tensor perturbation mode on the constraints of dark energy EoS. In the lower left 
panel of figure [U we present the constraint on the EoS parameters of dark energy from the models including tensor 
perturbation mode (green dotted lines). The tensor perturbation mode mainly contributes the CMB temperature 
anisotropics on the very large scales. Meanwhile, the dark energy component also affects the low multipoles temper- 
ature power spectrum of CMB through the late-time integrated Sachs- Wolfe (ISW) effect [sj- When including the 
tensor perturbation mode, we obtain the la constraints: wq = — 0.97± 0.18 and Wa = — O.ll^Q ^g. The constraints 
do not change significantly, considering the large error bars. 

Finally, we investigate the degeneracy between the EoS of dark energy and the massive neutrino. In the lower right 
panel of figure [U we show the constraint on wq and Wa when considering the massive neutrinos (blue dash-dot lines) . 
In this case, the constraints on the EoS parameters of dark energy are obviously weakened and the median values 
are shifted: wq = —0.92 ± 0.20 and Wa — —0.86 ± 1.00 (68% C.L.). This results in a strong anti-correlation between 
^ rrii, and the EoS parameters of dark energy, which mainly comes from the geometric feature of our universe fl(| . 
In addition, dynamical dark energy will modify the time evolving potential wells which affect CMB power spectra 
through the late time ISW effect. Dynamical dark energy can leave imprints on CMB, LSS power spectra, and Hubble 
diagram, nonetheless these features can be mimicked by cosmic neutrino to some extent [38|. In our analysis, we find 
that the correlation coefficients are 0.08 and —0.37 between '^rrii, and wq or Wa, respectively. 



2. Inflationary parameters 

Currently, the cosmological observational data are in good agreement with a Gaussian, adiabatic, and scale- invariant 
primordial spectrum, which are consistent with single-field slow-roll infiation predictions. Measuring the spectral index 
Us of primordial power spectrum is important for the studies of infiation paradigm. In the literature, the constraint 
on Ug is usually obtained by fitting the basic 6-parameter ACDM model. However, the constraints can be changed by 
the correlations between n^. and other cosmological parameters, such as the running of spectral index as, defined as 
as = dus/dhik, the tensor perturbation mode, as well as the EoS parameters of dark energy. 

To investigate these correlations, we constrain Ug in several different input cosmological models: the standard ACDM 
model which assumes a pure scale invariant power-law power spectrum (noted as ACDM), the standard ACDM model 
with running spectral index (noted as "ACDMSag"), the standard ACDM model with tensor perturbation mode 
(noted as "ACDM©r"), the standard ACDM model with running spectral index and tensor perturbation mode 
(noted as "ACDM © as © r" ) , the standard ACDM model with running spectral index and massive neutrinos (noted 
as "ACDM © as © m,y" ) , the standard ACDM model with tensor perturbation mode and massive neutrinos (noted 
as "ACDM©r © TOj/"). Furthermore, we have studied the correlations between the dynamical dark energy and 
inflationary parameters (noted as "CPL", "CPL © ag" and ^^CPL © r" models, respectively). We list the numerical 
results in table II and plot the 1-dimensional distributions of Us when using different input cosmological models. 

In the standard 6-parametcr ACDM model, we obtain the la constraint on Us is Us — 0.969 ± 0.011, which is 
consistent with previous works 2] . When considering the running of spectral index, the constraint of ris is significantly 
weakened, Us = 0.951 ± 0.020 (68% C.L.), while the 68% constraint of as is ^ -0.018 ± 0.016. The change of 
the constraint on rig is about 1.6 a. We also obtain that the correlation coefficient between Ug and as is 0.83, which 
implies the strong correlation between them. When including tensor perturbation mode, the constraint on Ug slightly 
changes: Ug = 0.974 ± 0.012 at the 68% confidence level. When varying both ag and r simultaneously, the 1 a error 
bars of Inflationary parameters are sig nificantly enlarged: Us = 0.946tJ^ !]^? and as = -0.038 ± 0.022. The correlation 
coefficients among ris, as and r are 0.70, 0.35 and —0.60, respectively. In figure[3]we plot the 2-dimensional constraints 
on the panel {ng,as) in the upper two panels with or without considering the tensor perturbation, respectively. 

The correlation between the early universe inflation and late time accelerating expansion can also affect the con- 
straints on inffationary parameters. In time-evolving dark energy model [CPL), we obtain the 68% C.L. constraint 
Ug = 0.966 ± 0.013, which is similar with that in the ACDM framework. However, when considering the running of 
spectral index, the constraints of Ug and as are obviously weakened, namely, rig = 0.9301^026' ^^'^ = ~0.029tQ 020 
at the 68% confidence level. More importantly, their mean values are shifted significantly, when comparing them from 
the standard ACDM model. In the "CPL © r" model, we obtain the constraints rig = 0.979to;oi8 (1 a) and r < 0.20 
(2 a). Their error bars become larger which is clearly shown in the lower right panel of figure [31 due to the correlation 
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TABLE II. la constraints on the Inflationary parameters Us, as, and r from Union2.1+WMAP7+BAO+IIST. For the 
weakly constrained parameters, we quote the 95% upper limits instead. 
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FIG. 2: 1-dimensional probability distribution of the spectral index Us in different input cosmological models. The black 
solid line is given by the standard ACDM model, while the red dotted, blue dash-dotted line, green dashed are obtained from 
the ACDM©Q!s, ACDM®r and ACDM®m^ models, respectively. For comparison, we also show the distribution from the 
standard CPL model (purple dash-dot-dot line). 



between the dark energy and the tensor perturbation we discuss before. 

We also study the correlations between the massive neutrino and rig, as, r. With the massive neutrino, we obtain 
the constraints Us = 0.960lg;Q23 i^'^) ^'^^ < ^^-^^ (^f) i^i the "ACDM©m.i^ © r" model, whose error bars are much 
larger than those in the massless neutrino case. 



3. Neutrino mass 



Neutrino, as a hot dark matter candidate, can be constrained from the free streaming modification of the transfer 
function of the matter power spectrum. Thus, the constraints on neutrino mass are correlated with the shape of 

TABLE III. The 95% confidence level upper hmits on the total mass of neutrinos from Union2.1+WMAP7±BAO±HST. 
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LCDM © 


< 0.45 eV 




LCDM ®m^®r 


< 0.55 eV 


-1.1 
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< 0.42 eV 


-1.1 
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FIG. 3: 2-dimensional cross correlation between tIs and Os in different input cosmological models: The black solid lines are given 
by the standard ACDM model, the red dotted lines, blue dash-dot lines and green dashed lines are given by ACDM (Bets (Br, 
ACDM © nil, ffi Us and CPL © as models, respectively. 




FIG. 4: 1-dimensional probability distributions of the neutrino mass in different input cosmological models. The black solid 
line is given by the standard ACDM model, while the blue dotted, red dashed and the green dash-dot lines are given by the 
ACDM ®m„ ® as, ACDM © m-i, © r and CPL © models, respectively. 



primordial power spectrum. Since the free streaming effects can be compensated by the modification of the primordial 
spectrum, the constraints on total neutrino mass can be influenced by the inflationary parameters, as we discuss above. 
Moreover, the behaviour of dark energy EoS can also affect the constraint of the total neutrino mass. In table III and 
figure|4]we show the numerical constraints on the massive neutrino from different input models: the standard ACDM 
model (noted as "ACDM©r7ii,"), the standard ACDM with tensor perturbation mode (noted as "ACDM® mi. © r"), 
the standard ACDM with running of spectral index (noted as "ACDM©to^ © and the CPL dynamical dark 

energy model (noted as ^^CPL © niu')- 

In the framework of standard ACDM model, the 2 a upper limit on the total neutrino mass is ^ < 0.45 eV. 
When including the tensor perturbation mode, the constraint is slightly weakened, namely, ^ TOi, < 0.55 eV at the 
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TABLE IV. The 1 a constraints on parameters from the simulated mock data with the fiducial value fifc = —0.013. 
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Wa 
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-0.0058 
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95% confidence level. As we know, both the negative as and the massive neutrino lead to a damped power on small 
scales. The effect of massive neutrinos may be compensated by a non- vanishing running of the primordial spectrum. 
On the other hand, the running is negative, this will lead to even more stringent constraints on the neutrino mass 
compared with fittings in the constant scalar spectral index cosmology jsoj . In practice, we find the 95% upper limit 
on the total neutrino mass becomes smaller, when varying the running of the spectral index, ^ nii, < 0.42 eV. 

The massive neutrino is strongly correlated with the EoS of dark energy. Therefore, the behaviour of dark energy 
model, especially the dynamical dark energy model, will significantly enlarge the constraint of total neutrino mass. 
In the "CPL © rriiy" model, we obtain the 2a upper limit on the massive neutrino is "Y^nii, < 0.81 eV, as shown in 
figure m 

B. Future Data 

From the results presented above, we see that there are some strong degeneracies among cosmological parameters, 
such as w and ^ m„, Ug and ag- Due to the precision of current observations, these degeneracies can not be broken 
and will weaken the constraints of cosmological parameters. Therefore, it is worthwhile discussing whether future 
observations could give more stringent constraints on the cosmological parameters and break these degeneracies 
efficiently. For this purpose, we have performed a further analysis and we have chosen the fiducial model in perfect 
agreement with current data: ^th'^ = 0.0227, fl^h'^ = 0.114, 1^ = 1.041, t = 0.0865, wq = -1, Wa = 0, Us = 1, 
As = 2.85 X 10-9 at K = 0.05 Mpc'^ 

1. Qk and w{z) 

We use the fiducial mock data to study the degeneracy between the EoS of dark energy and ftk- When we simulate 
the future data, we choose the fiducial value of curvature flk = —0.013 which is allowed by current observational data. 
We summarize the numerical results in table IV. 

Firstly, we use the correct input cosmological model "CPL © ilk" to constrain parameters. By combining the 
future mock data from PLANCK, Euclid and supernovae, we obtain the constraints on the EoS parameters: wq = 
— 1.0001qq280 ^^"^ ~ '^■00 i 0.130 at the 68% confidence level. When comparing the constraints from the current 
observations, the future data constrain the EoS parameters much more stringently. The error bars of wq and Wa are 
shrunk by a factor of 6. In figure [5] we plot the 2-dimensional constraint on the panel {wo,Wa) (black solid lines). 
Meanwhile, the mock data also give very tight constraint on the curvature: fifc = — 0.0130j;Q oo2o (68% C.L.). 

Besides the correct input model, we also use a biased model to fit the mock data, in order to study the degeneracy 
between ilk and w. Here, we assume the flat universe (ilk = 0). Due to the strong degeneracy between them, in 
this case the constraints of wq and Wa are quite different from those of the correct input model. The 1 a constraints 
are wq — ~l-068lQ g22o = 0.481^° q72- One can see that the error bars are smaller than above, due to less 

free parameters in this case. More importantly, the mean values of wq and Wa are shifted obviously, as shown in the 
red dashed lines of figure [SJ When considering the error bars, these mean values will be ruled out at more than 3 a 
confidence level, which implies that the future observations could break this degeneracy very well. 

We also use the constant w model {wa = 0) to fit the mock data again and plot the results in figure IH When 
varying flk, we obtain the constraints: w = -1.0000 ± 0.0085 and Qk = -0.0127io °oio (68% C.L.), which recover the 
fiducial model very well. However, when assuming the flat universe, due to the degeneracy, the 1 cr constraint on w 
is shifted: w = — 0.9407lQ QQ5g, which will be excluded by the future data at more than 6 a confldence level. In this 
case, the minimal are very large, which implies that this kind of model does not fit the mock data well. 
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FIG. 5: 1 and 2-dimensional constraints on the panel (ioo,u)a) from the simulated mock data with the fiducial value f2fc = — O.Of 3. 
The black solid lines are obtained by using the correct model which allows varying. For comparison, we also use a biased 
model which fixes the curvature Jlfc = to constrain the EoS of dark energy, shown as the red dashed lines. 
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FIG. 6: i-dimensional constraint on w and 2-dimensional constraint on the panel (VLh^vS) from the simulated mock data with 
the fiducial value Q.^ = — O.OfS. The black solid lines are obtained by using the correct model which allows ilfc varying, while 
the red star denotes the best fit model when assuming the flat universe. 



2. ^ rriv and w 

Similarly, we use the same method to study the degeneracy between w and ^ m^. When generating the simulated 
mock data, we choose the fiducial value of total neutrino mass J2 "^i' — 0.4 eV, which is in agreement with the current 
constraint at 2 cr confidence level. 

In figure [71 we plot the 2-dimensional constraints on wq and Wa obtained by using different input cosmological 
models. Firstly, we consider the correct model {"CPL © TOi/") and obtain very tight constraints on the EoS of dark 
energy and the total neutrino mass: wq = -1.000 ± 0.0240, Wa = 0.00 ±0.110 and Y.'^'y = 0-400to;°i°eV at the 68% 
confidence level. The future measurements could give conclusive conclusion about the total neutrino mass. However, 
if assuming the massless neutrino (^to^ = 0), wc will obtain the biased constraints on the EoS of dark energy: 
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TABLE V. The la constraints on parameters from the simulated mock data with the fiducial value = 0.4 eV. 



models 


Wo 




J2 


Ax' 


fiducial model 


-1 





0.4 eV 




CPL ® Tn„ 


-1.000 ±0.0240 


0.00 ±0.110 


400+° °^° 





CPL 


-1.031 ±0.0140 


0.589 ± 0.043 


set to 


40 


Constant w © 


l.UUi„g 0120 


set to 


0.400tSJ;?,i^ 





Constant w 


U.OOUO_o.0038 


set to 


set to 


104 



1 

0.8 

0.4 
0.2 



-1.05 -1 -0.95 

0.8 
0.6 
0.4 
g" 0.2 


-0.2 

-0.4 

-1.05 -1 -0.95 -0.4 -0.2 0.2 0.4 0.6 

w„ w 

a 

FIG. 7: 1 and 2-dimensional constraints on the panel (wo,Wa) from the simulated mock data with the fiducial value ^ = 
0.4 eV. The black solid lines are obtained by using the correct model which allows ^ varying. For comparison, we also use 
a biased model with the massless neutrino to constrain the EoS of dark energy, shown as the red dashed lines. 



Wo = —1.031 ± 0.0140 and Wa = 0.589 ± 0.043 (68% C.L.), due to the strong degeneracy between w and J2 The 
mean values arc significantly different from the input fiducial model, especially for the parameter Wa- Consequently, 
the minimal of this case will be very large, which implies that this best fit model has been ruled out by the data 
at more than 3 a confidence level. We list the numerical results in table V. 

Again, we also use the constant EoS of dark energy to study this degeneracy. In the framework of the correct 
w © nil, model, the la constraints are w — — l-00llgQ^2o and "^f = 0-400to'oi5- When forcing '^nii, = 0, the 
best fit model of w is shifted Xo w — — 0.8503lg'QQ3g, which is more than 10 cr away from the fiducial value of w, as 
shown in figure IH The future observations could break this degeneracy efficiently. 

3. Us, a. and r 

In order to study the degeneracies between Us and a^, r, we simulate two mock data sets: one with as = —0.05 
and the other one with r — 0.2, which are in agreement with the current data. In table VI we list the constraints of 
these Infiationary parameters from two mock data. 

Using the correct input model, we obtain very stringent constraints on Us and as', — l.OOOOiO.OOlO, as — 
— 0.0502lg;Qg4g (68% C.L.), whose error bars are almost shrunk by a factor of 10. When assuming the spectral index 
is scale-independent {as = 0), the constraint of Us changes: Us = 1.0052lg Qg}g at the 68% confidence level. The 
mean value of Us is shifted significantly, when comparing with the obtained error bar. We compare these two results 
in figurc[51 The obtained correlation coefficient between and as is 0.1, which means this degeneracy does not show 
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FIG. 8: 1-dimensional constraint on w and 2-dimensional constraint on the panel {w,^ m^) from the simulated mock data with 
the fiducial value = 0.4 eV. The black solid lines are obtained by using the correct model which allows varying, 

while the red star denotes the best fit model when assuming the massless neutrino model. 

TABLE VI. 1 a constraints on Inflationary parameters from the simulated mock data. 



models 


ris 


as 




fiducial model 


1 


-0.05 




correct model 


1.0000 ± 0.0010 


-0 n5n2+°°°^'-* 

U.UOU/_o 0040 





biased model 


1 00^2+°°°^^ 


set to 


78 


models 


Tis 


r 


Ax' 


fiducial model 


1 


0.2 




correct model 


1 0000+°°°" 

i. UUUU_Q 0010 


n 200+° °^^ 

U.ZUU_o.o20 





biased model 


0.9975±0.0009 


set to 


60 



in the future data significantly. 

The similar situation is also found by the study about the degeneracy between and r. In figure [TOl we show the 
constraints on and r from the correct and biased input models, respectively. When including the tensor perturbation 
mode, the la constraints are: Ug = l.OOOOlg'ooiJj!) and r = 0.200lg'Q2o- Neglecting the tensor perturbation, the 
constraint on Ug is shifted to a lower value: Ug = 0.9975±0.0009. We obtain the correlation coefficient between Ug 
and r is —0.05, which implies that the degeneracies between Ug and r is not notable for the future data sets. 

C. Conclusion and discussion 

In this paper we perform global fitting analyses on serval cosmological models and present the constraints on 
cosmological parameters from the latest astronomical observations, including WMAP7, BAO measurement from 
the SDSS-II survey and supernovae "Union2.1" sample. We find that, with different input theoretical models, the 
constraints on some cosmological parameters can be very different, such as the equation of state of dark energy, 
the total neutrino mass and the index of primordial power spectrum, due to the strong degeneracies among these 
parameters. 

With the high quality mock data, such as future SN, CMB and LSS observations, we further study the impact 
of these degeneracies on the constraints of cosmological parameters. Using these accurate future data, the obtained 
results show that, some degeneracies can not be neglected, such as the EoS of dark energy and the curvature, the EoS 
and the massive neutrino. Ignoring them forcibly will lead to seriously biased constraints. 
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FIG. 9: 1-dimensional constraint on Us and 2-dimensional constraint on the panel {ns,as) from the simulated mock data with 
the fiducial value a, — —0.05. The black solid lines are obtained by using the correct model which allows as varying, while the 
red star denotes the best fit model when fixing — 0. 
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FIG. 10: 1-dimensional constraint on Us and 2-dimensional constraint on the panel {ns,r) from the simulated mock data with 
the fiducial value r = 0.2. The black solid lines are obtained by using the correct model which allows r varying, while the red 
star denotes the best fit model when neglecting the tensor perturbation mode. 
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